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Robust negative differential conductance and 
enhanced shot noise in transport through a 
molecular transistor with vibration assistance 

Bing Dong, X. L. Lei, and N. J. M. Horing, 



Abstract — In this paper, we analyze vibration-assisted sequen- 
tial tunneling (including current-voltage characteristics and zero- 
frequency shot noise) through a molecular quantum dot with 
two electronic orbitals asymmetrically coupled to the internal 
vibration. We employ rate equations for the case of equilibrated 
phonons, and strong Coulomb blockade. We find that a system 
with a strongly phonon-coupled ground state orbital and weakly 
phonon-coupled excited state orbital exhibits strong negative 
differential conductance; and it also shows super-Poissonian 
current noise. We discuss in detail the reasons and conditions 
for the appearance of negative differential conductance. 



I. Introduction 

IN a nanoelectromechanical system, mechanical motion can 
affect electrical transport and vice versa. [1] This has 
stimulated a large body of experimental [2]-[8] and theoretical 
[9]— [19] work concerning phonon-assisted resonant tunneling 
through a quantum dot (QD) with strong coupling to an 
internal vibrational (phonon) mode (IVM), since tunneling of 
a single electron may induce a displacement of the mobile 
mechanical structure (which can be used as a prototype of a 
nanoscale displacement sensor). On the other hand, the vibra- 
tion of an adsorbed molecule can induce a change of tunneling 
current, providing a mechanism for sensing its adsorption 
and its identification. Current progress in nanotechnology has 
facilitated the fabrication of single-electron tunneling devices 
using molecules, for instance, a single C 6 o(i4o) molecule 
and carbon nanotubes (CNTs). In particular, CNTs are ideal 
systems for designing electromechanical devices since they 
have a small diameter, a low mass, hence a strong electron- 
phonon coupling (EPC) constant, and more importantly a 
relatively high quality factor, up to 10 4 . 

Recent experimental studies of a long suspended CNT sam- 
ple connected to two electrodes has revealed the appearance 
of negative differential conductance (NDC) at the onset of 
each phonon-assisted current step. [8] Theoretically, several 
previous calculations have been carried out for a system of a 
single level coupled to the phonon mode, and have revealed 
that the combined effect of low vibrational dissipation, i.e. an 
unequilibrated phonon (hot phonon), and strongly asymmetric 
tunnel-couplings to the left and right electrodes is responsible 
for the appearance of NDC at an appropriate EPC strength. [9], 
[15], [16], [20] However, this NDC in the single-orbital system 
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is quite weak and can be destroyed easily by a dissipative 
environment [20], and even by applying an ac-bias voltage. 
[21] Nonetheless, it is desirable to search for a strong and 
robust NDC property in such transistors by tailoring material 
parameters for possible application in functional devices. 

Our recent work has examined unequilibrated vibration- 
assisted tunneling through a molecular QD with two electronic 
orbitals having asymmetric couplings to the electrodes, with 
both strongly interacting with an IVM. [20] We found a strong 
NDC and correspondingly enhanced current noise due to a 
trapping effect, and we analyzed in detail the influence of the 
uneqilibrated vibration on its appearance. On the other hand, 
Nowack and Wegewijs have considered vibration-assisted tun- 
neling through a two-level QD with asymmetric couplings to 
an IVM. [17] Their results display strong NDC through a 
competition between different Franck-Condon (FC) tailored 
tunneling processes of the two levels, which is completely 
robust even against strong relaxation of the vibration. In this 
paper, we further investigate tunneling in this two-orbital 
system with asymmetric couplings to the vibration in the 
equilibrated vibration condition, focusing on the reasons and 
conditions for the appearance of strong NDC, as well as its 
corresponding zero-frequency shot noise. We find that the 
current noise is always greatly enhanced (having a giant Fano 
factor) for this kind of system. 

The outline of this paper is as follows. In Sec. II, we 
describe the model system that we study. We also exhibit the 
rate equations in a number-resolved form at high temperature 
to describe phonon-assisted resonant tunneling in the limit of 
strong relaxation (equilibrated phonon). In Sec. Ill, we then 
investigate in detail the vibration-assisted transport and shot 
noise properties of a molecular two-orbital QD. Finally, a brief 
summary is given in Sec. IV. 

II. Model and rate equations 

Here we consider a generic model, as shown in Fig. 1, for 
a molecular QD with two spinless levels, one as the highest- 
occupied molecular orbital (HOMO), £i, and the other as 
the lowest-unoccupied molecular orbital (LUMO), £2, coupled 
to two electrodes, left (L) and right (R), and also linearly 
coupled to an IVM of the molecule having frequency uj n 
with respective coupling strengths Xi and A 2 . The model 
Hamiltonian is 

H = Hi eads + H mo i + H T , (la) 
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with 



H, 



leads 



E £ ^ C l 



77k 77k' 



(lb) 



77, k 



= e J c j c j + Unin 2 



.3=1,2 _^ 

+w a t a+ Aj-cjcj-^ + a), (lc) 

.7=1,2 

#t = ^(^c^+H.c), (Id) 

»7,kj 

where c^ k (c^k) is the creation (annihilation) operator of an 
electron with momentum k, and energy e^u in lead 77 (77 = 
L,R), and cj (c,) is the corresponding operator for a spinless 
electron in the jth level of the QD (j = 1, 2). [/ denotes 
interdot Coulomb interaction and n, 



have poo, P11, and P22 denoting probabilities for the electronic 
states: zero-electrons in the dot (0), single-electron occupation 
in the HOMO (1) or the LUMO (2), respectively. 

In the weak tunneling regime and high temperature approx- 
imation, T V j <C T, u>o (T V j is the tunneling rate between 
lead 77 and level j, and T is the temperature), we have 
established generic rate equations for these electron occupation 
probabilities including the role of phonon-assisted resonant 
tunneling in the eqilibrated phonon condition employing a 
quantum Langevin equation approach. [20] In order to employ 
MacDonald's formula for calculating shot noise [23], [24], we 
write these rate equations in a number-resolved form: 



•(0 
Poo 



r- o {l ~ 1] 



is the electron 



P11 
.(i) 
P22 



-(rf + 

r+ n {l) + r+ n {l+1) -r~o K ' 

1 LlPoO ^ 1 RlPOQ 1 1 Pi 



+ ^ L2P2I 



• r- 

1 #2^22 



,+ (I) , p+ (1+1) 
L2P00 ' 1 R2P00 



r -J0 

1 2 P22 1 



(3a) 
(3b) 
(3c) 



where (t) represents the electron occupation probability 



number operator in level j. In this paper, we consider strong 
Coulomb repulsion, U — > 00, which suppresses simultaneous 
occupation of the two orbitals. (a) is phonon creation 

(annihilation) operator for the IVM. Xj represents the EPC at time i" together with I electrons arriving at the right lead 
constant of an electron in level j; V nj describes the tunnel- due to tunneling events. Obviously, Pjj (t) = T,iPjj( t )- The 
coupling between electron level j and lead 77. We use units 
with h = ks = e = 1 throughout the paper. In addition, 
it is worth noting that we do not consider the EPC-induced 
coupling between the two MO's in the Hamiltonian, Eq. ( fTc| , 
because it is much weaker than the coupling of a single MO. 



electronic tunneling rates are defined as 



r , + = E r ™-=E r wE^(^+nwo). 



>IJ 




r- = Vr 

n 

= E r '"' E £ « t 1 ~ A ( ? i ~ ULU °^ ' 



(4a) 



(4b) 



T) n 



Ci = e-^ 2n ^e n ^ T l n (2g^n B (n B + l)),(4c) 

where /,,(e) = [1 + e( e ~ l * r "' T ]~ 1 is the Fermi-distribution 
function of lead 77 {/j, n is the chemical potential of lead 77), 



n B 



Fig. 1. Schematic diagram of the model system. 



It is well-known that the electron-phonon interaction term in 



is the Bose distribution function and 
is the nth Bessel function of complex argument. 
The current flowing through the system can be evaluated by 
the time rate of change of electron number in the right lead 

as 

,1 , — , 

(5) 



Eq. ( lc 1 can be eliminated by a canonical transformation, [22] 



I=Nn(t) = ±J2 lp(l) W 
1 



H = e^ He~ s , with S — (gini+g2n 2 )(a^ — a) (c/j — Xj/wq), 
leading to a renormalization of the energy parameters, ej — 
Sj — 9jXj, and of the tunnel-coupling, V„j — V„j exp[<?j(a^ + 
a)]. The transformed Hamiltonian is 



where 



(I), 



,(0/ 



(6) 



Ht — 



Efai- 



V/k C 



H.c), 



(2a) 



(2b) 



is the total probability of transferring I electrons into the right 
lead by time t. The zero-frequency shot noise with respect to 
the right lead is similarly defined in terms of P"'(t) as well: 
[20], [23], [24] 



S(0) 



> Ai A2 



with U=U -2- 

In this paper, we consider the case of strong relaxation 
of the IVM. In this situation, the vibration is always in an 
equilibrated state during tunneling processes, i.e., the excited 
phonon relaxes very quickly due to strong dissipation, before 
the next electronic tunneling event takes place. As a result, 
in this strongly Coulomb-interacting two-level QD system, we 



dt 



(7) 



To evaluate S(Q), we define an auxiliary function Gjj(t) as 



= EW) 



(8) 



whose equations of motion can be readily deduced employing 



the number-resolved QREs, Eqs (3ai-(3ci, in matrix form 
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G(t) = MG{t)+Qp{t) with G(t) = (G 00 ,Gii,G 22 ) T and 
p(t) — (poo, Pn, P22) T ■ Applying the Laplace transform to 
these equations yields 



G(s) = {sI-M)- l Gp{s), 



(9) 



where p(s) is readily obtained by applying the Laplace 
transform to its equations of motion with the initial condition 
p(0) = p s t (Pst denotes the stationary solution of the rate 
equations). Due to the inherent long-time stability of the 
physical system under consideration, all real parts of nonzero 
poles of p(s) and G v (s) are negative definite. Consequently, 
the divergent terms arising in the partial fraction expansions 
of p(s) and G' n (s) as s — > entirely determine the large- 
t behavior of the auxiliary functions, i.e. the zero-frequency 
shot noise, Eq. Q. This scheme facilitates evaluation of the 
zero-frequency shot noise as a function of bias-voltage for 
equilibrated-vibration-assisted tunneling through a molecular 
QD. 

III. Results and discussions 

In this section, we proceed with numerical calculations of 
the current /, the zero-frequency Fano factor F = S(0)/2I for 
the two-orbital model. We set the symmetric tunnel-couplings 
as T L1 = T m = T L2 = T R2 = O.OIwq, take a weak EPC 
constant of the excited level, g 2 — 0.1, with the energy gap 
between the LUMO and HOMO as A = e 2 - £i = w /2, 
in our calculations (Here, several values of g\ are employed). 
Throughout the paper, the energy of IVM, loq, is chosen as 
the energy unit, and we set the temperature as T = 0.05o;o 
unless otherwise indicated and assume that the bias voltage V 
is applied symmetrically, hl = —[Ir = V/2. We have checked 
our parameter choices by comparing with the experimental 
data. [2], [8] Our numerical calculations yield the maximum 
current through the molecule as I max ~ 0.5T (see below), 
and the experimental result is about I max ~ 100 pA, therefore 
we obtain T ~ 1 //eV, which is nearly ~ O.Ol^o since the 
characteristic energy loq of the oscillation is found to be ~ 
0.5 meV for the experimental setup of Sapmaz, et al. [8] 

Figure 2 exhibits the calculated results. For comparison, 
we also plot the corresponding results of the single-level QD 
(which has only the ground orbital involved in transport) 
in Fig. 2(a,b). It is clearly evident that for the single-level 
system, the strong coupling between electronic and vibrational 
degrees of freedom induces a significant current suppression 
at low bias voltage due to a phonon-modified tunneling rate, 
and the noise is always sub-Poissonian, albeit the current is 
also largely suppressed by the strong EPC. Interestingly, it is 
seen that the current noise has an oscillatory feature with an 
approximate period of V = 2.0cc>o, if £i ^ ±nwo/2 (n is an 
integer), which can be explained based on Eqs. ([3a]>-(|4c]>. For 
the sake of simplicity, we assume low temperature, T — 0, in 
the following analysis. 

To start, we consider a system with |ei| < Wo/2, for 
example, E\ — — 0.25wo- Within the bias voltage window 
0.5u>o < V < l.5u! , we have the tunneling-in rate, x' = 
= TliC^, and the tunneling-out rate, x = = TmC^ 
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Fig. 2. Current I (a,c) and Fano factor F = 5(0) /2I (b,d) vs. the bias 
voltage for a single-level system, ei/u>o = —0.25 (a,b), and for a two-level 
system, ei/u>o = —0.25, A/tJrj = 1/2, gi = 0.1 (c,d), with various values 
of g\. Inset: Fano factor vs. bias voltage for the single-level system with 
gi = 2.5 and various energies of the QD level, e \ = (solid line), —0.25 
(dashed line), —0.5 (dotted line), and —0.75 (dotted-dashed line). 



with symmetric tunnel-couplings under large bias voltage. Its 
typical Fano factor is F\ = 1/2. However, when the bias 
voltage becomes 1.5wo < V < 2.5u>q, the n = 1 channel 
is open, which leads to x' = = + C\) but still 

has x — TZ — TjiiCq. The equivalent model changes to an 
asymmetric single-level system with a Fano factor given by 

x' 2 -I- x 

(10) 



(x' + x)- 



< 1. 



If the bias voltage increases further to 2.5wo < V < 3.5ujq, 
T]j becomes equal to again and the Fano factor changes 
back to 1/2. Therefore, because of the specific forms of 



(r 



0). This is equivalent to a single-level QD 



the phonon-assisted tunneling rates, Eqs. (4b I and (4c i, the 
single-level QD varies between two equivalent models, the 
symmetric model and asymmetric one, with increasing bias 
voltage, which generates the oscillatory feature of the Fano 
factor. Similarly, it is easy to verify that a system with e\ = 
always corresponds to the symmetric model with a constant 
Fano factor, Fi = 1/2. 

A system with |e"i| > luq/2 always corresponds to the 
asymmetric model, whose noise exhibits step structure only 
if |ei| = nuio/2 (n is a positive integer), but it is oscillatory 
otherwise, as shown in the inset of Fig. 2(b). Here, for 
illustration, we take E\ — — 0.5wo an d — 0.75wo as examples. 
Their zero-temperature tunneling rates are listed in Table I 
for three different bias voltages. It is obvious that the rates 
of tunneling-in and -out, x' and x, of the system with e\ — 
— 0.5wo experience a simultaneous increase with increasing 
bias voltage, leading to a step-down of the current noise; 
but their increases become out of phase with respect to bias 
voltage for a system with e\ = — 0.75wo^ the tunneling-in 
rate (x') suffers an increase first and the tunneling-out rate 
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(x) remains unchanged from V = 2. Own. to 3. Own, then x 
increases at V = 4. Own, which generates an oscillation of the 



noise [Eq. (10 1]. 

In the case of a single-level QD, there is no NDC and no 
super-Poissonian noise. However, the situation is different for 
two-level systems. We find a strong NDC and a huge Fano 
factor (F ~ 100 for gi = 2.5) with e x = -0.25w , A = w /2, 
and gi > 1.0, as shown in Figs. 2(c,d). Actually, we find from 
Fig. 3 that no NDC occurs only when e\ — ±nu>o/2 between 
— 1.5wq < £ i < l.Own. 




Fig. 3. The differential conductance, dl j dV (ei, V) for gi = 2.5, <?2 = 0.1, 
and A = lu /2 at T = 0.05a;o. 



This peculiar feature can be understood from the bias- 
voltage-dependent electron occupation probabilities. For a 
system with e\ = — 0.25wn [Fig. 4(a)], considering increasing 
bias voltage up to 0.5wn, the excited orbital (LUMO) starts 
to be occupied and then contributes to current. When the 
bias voltage arrives at 1.0a>o> only two states, In and 2o (i n 
denotes the n-phonon-assisted tunneling channel of orbital 
i), are involved at the bias voltage window with respective 
symmetric tunneling rates, x = Pj^ = = T L1 C\ <C T L1 
and y = T+ 2 = T R2 = T L1 C 2 sa T L1 (due to g 2 < 
1 -C gf). This situation is equivalent to a two-level QD with 
symmetric tunnel-couplings under a large bias voltage at zero 
temperature, as indicated in Fig. 5(a). Thus, we have attained 
the maximum value of pu (pn ~ p no ~ P22 — 1/3) and the 
current and Fano factor are given by 

x + Vi 



-Tli, 

2(x 2 + y 2 ) + xy 
9xy 



(11) 



(12) 



respectively. From Eq. ( [T2| , we obtain the huge Fano factor 
(due to 1 « y), F a ~ 2y/9x ~ 114. If the bias voltage 
increases further, for example to 1.5uiq, the tunneling channel 
li is in alignment with the Fermi energy of the left lead. 
As a result, an electron in the QD prefers to re-occupy the 
ground orbital, which causes a suppression of current due 
to its stronger EPC constant than that of the excited orbital, 
9i ^ #2- In this case, the zero-temperature equivalent model 
is depicted in Fig. 5(b), where tunneling through the ground 
orbital becomes asymmetric, x' = = Tli(£q + C\) > 



x = T m . Solving the simplified rate equations, one obtains 

x' x 



P11 
h 



x' + 2x 
x(x' + y) 
x' + 2x 
2x'(y 2 + 



> P22 = Poo = 



t 2 — xy) + y(x' 2 



2x' 



2x 2 ) 



y(x' + 2xf 



(13) 
(14) 

(15) 



It should be noted that if the bias voltage increases further 
to 2.5wo < V < 3.5wo, the system changes back to the 
symmetric model (a) with a stronger tunneling rate, x = 
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Fig. 4. (a,b) Electron occupation probabilities, (c) current /, and (d) Fano 
factor F vs. the bias voltage for the two-level system with A/i^o = 0.5 and 
gi = 2.5. (a) is plotted for i?i/ojo = —0.25; (b-c) are for £i/aj = —0.5 
(thick lines), (thin red lines), and —0.75 (thin blue lines). 




(b) 




x' 



Fig. 5. Schematic diagram of two equivalent models for a two-orbital system. 
The arrows denote the directions of electron tunneling. 

Similar to the single-orbital QD, if the energy of the ground 
orbital is zero, e\ — 0, the system always corresponds 
to the equivalent model (a). As a result, it exhibits step 
structures in the electron occupation probabilities [Fig. 4(b)] 
and current [Fig. 4(c)]. However, a system with |ei| > wo/2 is 
always equivalent to model (b): (1) if |ejJ = nujo/2, electron 
occupation probabilities also exhibit step structures because of 
simultaneous increases of x' and x (Table I) with increasing 
bias voltage. Correspondingly, no NDC occurs in its current; 
(2) On the contrary, if 7^ nuio/2, the tunneling-in rate 
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TABLE I 

Tunneling rates, x' = r£j and x = T^, of the ground orbital with ex = -0.5w and -0.75oj at V/lu = 2.0, 3.0, and 4.0, 

RESPECTIVELY, AT ZERO-TEMPERATURE. 
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x' increases first with increasing bias voltage, which means 
that an electron has a greater opportunity to enter the ground 
orbital, i.e., p\\ increases. With further bias voltage increase, 
the tunneling-out rate x also increases, leading to a decrease of 
Pix- Therefore, this system also exhibits oscillatory behavior 
in its electron occupation probabilities, and thus exhibits an 
NDC feature in its current. The noise in all these cases is still 
super-Poissonian due to x' .x ^ y [Fig. 4(d)]. 
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Fig. 6. Current I (a,c) and Fano factor F (b,d) vs. bias voltage for the 
two-level system with ei/oJo = —0.25 (a,b) and Ex/^o = —0.75 (c,d) 
for various temperatures T. Other parameters: gx = 2.5, g2 = 0.1, and 
A/wo = 0.5. 

Finally, we exhibit the temperature dependence of the NDC 
and super-Poissonian current noise for the two-level systems 
in Fig. 6. It is found that the NDC behavior of the current 
gradually becomes weaker with increasing temperature and 
nearly disappears at a relatively high temperature, T = 0.5wo, 
but the noise still remains strongly enhanced, F ^> 1, in 
the low bias-voltage region even though the Fano factor F 
decreases with increasing temperature. 

IV. Conclusions 

In summary, we have analyzed equilibrated phonon-assisted 
tunneling and its zero-frequency current noise through a single 
molecular QD having two electronic orbitals with asymmetric 
couplings to the IVM using generic rate equations. Our results 
show that due to asymmetric EPC, the two-orbital molecular 
transistors exhibit strong and robust NDC and enhanced shot 



noise with a large Fano factor. We have also fully analyzed 
the reasons and conditions for the appearance of NDC in these 
systems. 

Acknowledgment 

This work was supported by Projects of the National 
Science Foundation of China, the Shanghai Municipal Com- 
mission of Science and Technology, the Shanghai Pujiang 
Program, and Program for New Century Excellent Talents in 
University (NCET). 

References 

[1] M.L. Roukes, Phys. World 14, 25 (2001); L.Y. Gorelik, A. Isacsson, M.V. 
Voinova, B. Kasemo, R.I. Shekhter, and M. Jonson, Phys. Rev. Lett. 80, 
4526 (1998); A. Erbe, C. Weiss, W. Zwerger, and R.H. Blick, Phys. Rev. 
Lett. 87, 096106 (2001). 
[2] H. Park, J. Park, A. Lim, E. Anderson, A. Allvisatos, and P. McEuen, 

Nature 407, 57 (2000). 
[3] J. Park, A.N. Pasupathy, J.I. Goldsmith, C. Chang, Y. Yaish, J.R. Petta, M. 
Rinkoski, J.R Sethna, H. Abruna, PL. McEuen, and D.C. Ralph, Nature 
417, 722 (2002). 



and D. Natelson, Phys. Rev. Lett. 93, 266802 (2004); L.H. Yu and D. 
Natelson, Nano Lett. 4, 79 (2004). 

A. N. Pasupathy, J. Park, C. Chang, A.V. Soldatov, S. Lebedkin, R.C. 
Bialczak, J.E. Grose, L.A.K. Donev, J.R Sethna, D.C. Ralph, and PL. 
McEuen, Nano Lett. 5, 203 (2005). 

B. J. LeRoy, S.G. Lemay, J. Kong, and C. Dekker, Nature 432, 371 (2004); 
B.J. LeRoy, J. Kong, V.K. Pahilwani, C. Dekker, and S.G. Lemay, Phys. 
Rev. B 72, 75413 (2005). 

[8] S. Sapmaz, P. Jarillo-Herrero, Ya.M. Blanter, C. Dekker, and H.S.J, van 
derZant, Phys. Rev. Lett. 96, 26801 (2006); S. Sapmaz, P. Jarillo-Herrero, 
Ya.M. Blanter, and H.S.J, van der Zant, New J. Phys. 7, 243 (2005). 

[9] D. Bose and H. Schoeller, Europhys. Lett. 54, 668 (2001). 

[10] A.S. Alexandrov and A.M. Bratkovsky, Phys. Rev. B 67, 235312 (2003). 

[11] K.D. McCarthy, N. Prokof'ev, and M.T. Tuominen, Phys. Rev. B 67, 
245415 (2003). 

[12] A. Mitra, I. Aleiner, and A.J. Millis, Phys. Rev. B 69, 245302 (2004). 
[13] J. Koch and F. von Oppen, Phys. Rev. Lett. 94, 206804 (2005). 
[14] J. Koch, M.E. Raikh, and F. von Oppen, Phys. Rev. Lett. 95, 56801 
(2005). 

[15] J. Koch and F. von Oppen, Phys. Rev. B 72, 113308 (2005). 
[16] A. Zazunov, D. Feinberg, and T. Martin, P hys. Rev. B 7 3, 1 15405 (2006). 
[17] K.C. Nowack and M R. Wegewijs, cond-mat/0506552 (2005) 
[18] M R. Wegewijs, K.C. Nowack, New J. Phys. 7, 239 (2005). 
[19] F. Haupt, F. Cavaliere, R. Fazio, and M. Sassetti, Phys. Rev. B 74, 
205328 (2006). 

[20] X.Y Shen, Bing Dong, X.L. Lei, and N.J.M. Horing, Phys. Rev. B 76, 
115308 (2007). 

[21] Bing Dong, X.L. Lei, and N.J.M. Horing, Appl. Phys. Lett. 90, 242101 
(2007). 

[22] GD. Mahan, Many-Particle Physics. (Third edition, Kluwer Aca- 
demic/Plenum Publisher, New York, 2000). 
[23] D.K.C. MacDonald, Rep. Prog. Phys. 12, 56 (1948). 
[24] L.Y. Chen and C.S. Ting, Phys. Rev. B 46, 4714 (1992). 



